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Quantum Simulation of Spin Chains Coupled to 
Bosonic Modes with Superconducting Circuits 


U. Las Heras^ L. Garcia-Alvarez^ A. Mezzacapo^’^ E. Solano^’^ and L. Lamata^ 


Abstract We propose the implementation of a digital quantum simulation of spin 
chains coupled to bosonic field modes in superconducting circuits. Gates with high 
fidelities allows one to simulate a variety of Ising magnetic pairing interactions with 
transverse field, Tavis-Cummings interaction between spins and a bosonic mode, 
and a spin model with three-body terms. We analyze the feasibility of the implemen¬ 
tation in realistic circuit quantum electrodynamics setups, where the interactions are 
either realized via capacitive couplings or mediated by microwave resonators. 


1 Introduction 

A two-level system coupled with a single radiation mode is modeled by the ubiqui¬ 
tous and paradigmatic quantum Rabi model [HI that describes the most fundamen¬ 
tal interaction between quantum light and quantum matter. There have been many 
efforts, both in theory and experiments, to capture the features of this model in dif¬ 
ferent quantum technologies (2 El . These analysis have an impact on understanding 
about different quantum phenomena (H [5] [6l [7] [H [9l . 

The concept of a quantum simulator can be attributed to Feynman cni, and it 
refers to a controllable quantum platform that mimics the behaviour of another quan¬ 
tum system. Analog quantum simulators have been proposed and implemented in 
several quantum technologies, such as trapped ions HU [13, ultracold atoms C3, 
or superconducting circuits CdlElIIllEI- Digital methods based on discrete-time 
gate sequences ca in order to simulate dynamics of quantum systems have been 
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proposed and realized in trapped ions im, photonic systems 1^ . spin-photon hy¬ 
brid systems 1211 and superconducting circuits ll^l^[24ll^[^IZ7l . 

Here, we analyze the quantum simulation of arbitrary and generic models, where 
spin chains alone or coupled to bosonic modes are emulated in superconducting 
circuits 1^ . We use digital techniques in order to imitate systems whose dynamics 
in principle may differ from the ones of the experimental setups. Finally, we study 
the feasibility and efficiency of the implementation of three generic models in a 
realistic circuit quantum electrodynamics setup. 


2 Digital Quantum Simulations 


The goal of simulating diverse and generic models involving spin interactions and 
bosonic modes leads us to consider digital techniques, due to their suitability and 
flexibility for mimicking different dynamical structures. Hamiltonian dynamics can 
be approximated by the digital decomposition of the exact unitary evolution into 
discrete stepwise unitary operations, implemented by using quantum gates in an ef¬ 
ficient way HFll^. Digital methods are based on the Trotter formula, which allows 
us to expand the evolution operator of Hamiltonians that are written as a sum of 
terms, H = Y!j=\ Hj, into a product of evolution operators for the interactions given 
by the summands of the Hamiltonian, Hj. The Trotter expansion can be written as 


-iHt 


= ( .. . ^-iHNt/s\ ^ ^ 

i<j 


2s 


-'Em, 

k=3 


( 1 ) 


where the total time of the simulated dynamics is divided into 5' intervals in which 
the evolution associated to each summand of the complete Hamiltonian are applied. 
The error scales with t^/s for short times, as can be observed in the second term, and 
the upper bound for higher order error contributions is s\\Ht/s\\^y^^/kl > ||£’(k) ||sup. 

Our goal is to propose a systematic procedure using digital methods for simu¬ 
lating efficiently different models, namely spin-spin interaction and spins coupled 
to bosonic modes. First, we employ gates that commute with each other and do not 
produce digital error. For those that do not commute, we apply several Trotter steps 
because the more Trotter steps one applies, the smaller the digital error produced 
is. In realistic experiments, one has to take into account decoherence times and gate 
errors. Therefore, we have to regulate the number of steps in order to be able to 
perform the simulation before decoherence effects take place, and in order to reduce 
the accumulated gate error. Consequently, once the digital error is small enough ap¬ 
plying a certain number of Trotter steps, the error coming from the experimental 
setup always must be smaller than the digital one. 
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3 Quantum simulation of spin chains coupled to bosonic modes 
with superconducting circuits 


In this section, we present a method to implement the dynamics of several spin mod¬ 
els, coupled with bosonic modes, in circuit quantum electrodynamics setups. For 
this purpose, we take under consideration two different architectures of supercon¬ 
ducting circuits. We show how to use linear arrays of superconducting qubits with 
capacitive coupling between nearest neighbors 1^ to simulate the Ising model with 
transverse field. Then we simulate the behavior of a spin-chain coupled to a bosonic 
mode via a Tavis-Cummings interaction ED. Moreover, we show how to imple¬ 
ment many-body spin dynamics with highly nonlinear terms in superconducting 
qubits coupled to transmission line resonators acting as a quantum bus ||32l. 

In the following, we propose digital quantum simulations based on quantum gates 
implemented in superconducting architectures. Capacitive coupling setups allow 
one to implement ZZ gates, exp(—/0(7y (7|), for nearest-neighbor superconducting 
qubits by the sequence of two single qubit rotations along the z axis, Z(0), and a 
c-phase gate, CZ(0), as shown in Fig.[^ where 



/I 00 0 \ 

0 10 0 
00 1 0 
\0 0 0e-^^^ 


( 2 ) 


The current achievable fidelities in superconducting circuits are of 99.9% 
and 99.4% for the single and two-qubit (CZ) gates, respectively. They enable cir¬ 
cuit QED setups to be great candidates for digital quantum simulators where the 
stroboscopic application of gates is necessary. Notice that ZZi 2 { 6 ) = (Zi(0) C) 
Z2(0))CZi2(0) for 0 = 0/2. 


Fig. 1 Protocol for decom¬ 
posing a ZZ interaction 
between two qubits using 
single-qubit rotations, Z, and 
a c-phase gate, CZ. 



The use of quantum buses allows for the implementation of multi-qubit gates 
and spin-boson interactions, coupling the electromagnetic field in the resonator with 
superconducting circuits 1^1331 [Ml [35]| . 
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3.1 Ising model with transverse field via capacitive 
nearest-neighbour gates 

One of the most studied spin models in condensed matter is the Ising model with 
a transverse field, which is used for describing the behavior of interacting nearest- 
neighbor dipoles in the presence of a transverse magnetic field. The Hamiltonian of 
N spins can be written as 


Hitf = •/ L B L cTj-, (3) 

{jk) j 

where is the Pauli operator acting over the j-th spin with j = 1,...,A^, in the 
direction a = x^y^z. J stands for the coupling between nearest-neighbor spins and 
B is the coupling between a spin and the transverse field. Depending on the sign 
of J the model is ferromagnetic (/ < 0) or antiferromagnetic (/ > 0). In order to 
reproduce this interaction in superconducting circuits, we make use of a high-fidelity 
set of gates, as introduced in Eq. single-qubit rotations along the v direction, 
Xy(0) = exp(—/0(7y), and two-qubit ZZ gates, ZZjk{6) = exp(—/0cjy(7|). 

As shown in Section it is possible to decompose a complex interaction into 
discrete series of gates through a Trotter expansion. In order to implement the spin- 
spin interaction, we need to execute (A — 1) two-qubit gates. In this case, there is 
no digital error because all the gates in this decomposition commute, 

exp(-/ e Y, e-'e<T|(T| ... ^ 4 ^ 

iJk) 

with 6 =Jt,t being the simulation time of the experiment.The coupling among the 
spins and the transverse field can be simulated in a similar way using N single qubit 
rotations, 

exp(-/ (5) 

J 

with (j) =Bt. Given that the two interactions in Eqs. Q and ^ do not commute, one 
has to implement them in sequential short-time Trotter steps to minimize the digital 
error. In Eig.[^ we show a scheme of the protocol for the quantum simulation of the 
transverse field Ising model for four spins. The recent achievement of high-fidelity 
single and two-qubit (CZ) gates with superconducting circuits will allow one to 
perform many Trotter steps for several qubits, using hundreds of gates. 


Fig. 2 Protocol for digital 
quantum simulation of the 
Ising model with transverse 
magnetic field in terms of ZZ 
two-qubit gates and single 
qubit rotations along x axis. 
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In order to reduce the digital error, it is necessary to increase the number of Trot¬ 
ter steps. In Fig.|^ we plot a numerical simulation of the Ising model with transverse 
field for different digital steps. The simulated dynamics with digital decomposition 
is more accurate when compared with the exact dynamics when the number of Trot¬ 
ter steps is increased. 



Fig. 3 Digital simulation of the ferromagnetic Ising model with a transverse field for four spins 
in superconducting circuits, up to a phase of 0 = 4, with J/B = 2. The plot shows the fidelity 
of the digitally evolved state versus the ideally evolved one for different number of Trotter steps, 
^ = 6,8,10. The inset shows the overlap between the ideally evolved state with the initial state, that 
is, all qubits in |0)^. 


3.2 Spin chain coupled to a bosonic field mode via 
Tavis-Cummings interaction 

We now analyze a model consisting of a spin-chain with nearest-neighbour cou¬ 
plings interacting with a bosonic mode. In this sense, both free energies of the 
bosonic mode and spins are taken into account, as well as spin-spin and spin-boson 
interactions. The spin-spin evolution is modelled with the Ising dynamics, while the 
Tavis-Cummings model describes the interactions between spins and bosons. The 
resulting Hamiltonian is 

Hjtc = ft) ^ y (7? - y (7?g (acTt + «■'■ err) . 

j (jk) j 


(6) 
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Following the notation presented above, Oy is the Pauli operator along z direc¬ 
tion, is the creation(annihilation) spin excitation operator acting on the 

i-th spin and a{a'^) is the annihilation(creation) operator of the bosonic mode. CO 
and Q are the free energies of each boson and spin, respectively. Moreover, J is the 
coupling constant between nearest spins and g stands for the coupling among spins 
and bosonic field. 

The implementation in circuit QED requires the simulation not only of the spin 
dynamics, as in the previous example, but also of the bosonic mode. To achieve this, 
the superconducting qubits play the role of spins while the photons in a transmission 
line resonator emulate the bosonic excitations in the simulation. In order to perform 
the interactions of the model, it is necessary to couple the resonator to all the su¬ 
perconducting qubits. The Tavis-Cummings interaction appears straightforwardly in 
circuit QED setups once the rotating wave approximation is performed. 

Hi = (0i o'a + + gY^{aaf +a'(yj~), (7) 

j j 

being coi the frequency of the photons in the resonator, Qi the frequency associated 
with the excitation energy of the superconducting qubits, and g the qubit-resonator 
coupling constant. The spin-spin interaction for qubits j and k is achieved by means 
of the ZZ gate presented in Eq. Detuning to high frequencies the qubit-resonator 
interaction we are able to reproduce the model 

H{j,k) = (o' a^a + Y^ cr?-/crj. (8) 

j 

Since ^k')] = 0 V 7 ,/,^,^', we can define and implement sequentially 

the interaction 


H2= YH{i,k) = Oh a + Y^(^] - ( 9 ) 

{jk) j {jk) 

where (O 2 = {N — \)co' , 0.2 = and N the number of simulated spins, and 

it fulfills the condition exp(—fiE/ 2 ) = Il{jk) being t the execution 

time. 

Summing the interactions Hi and H 2 we recover the model we wanted to repro¬ 
duce, Hjtc- Nevertheless, [Hi^H 2 ] 7 ^ 0, so we need to employ the Trotter method in 
order to make the digital error decrease, as shown in Fig.[^ Moreover, for consider¬ 
ing the resonator photonic leakage, we have calculated the evolution of the system 
making use of the master equation, 

p = -i[Ht,p]-^KL{a)p, (10) 

where L{a)p = {2apa^ —a^ap — pa'^a) 12 is the Lindblad superoperator acting on 
a, K is the decay rate of the resonator, and Ht = {Hi^H 2 } is the Hamiltonian that 
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governs the evolution. Notice that we have considered a coherence time much longer 
for the qubits than for the resonator. In Fig. we plot the steps for implementing 
the protocol for four spins interacting with a bosonic mode. 



e 


Fig. 4 Fidelity F of the simulation of a four-spin chain coupled to a bosonic mode with circuit QED 
for different Trotter steps, ^ = 3,4,5. The upper curves correspond to larger number of Trotter steps. 
Here, the parameters of Hamiltonians Hi and H 2 are coi =271 x 200 MHz, = 27r x 180 MHz, 
g = 27r X 80 MHz, CO 2 = 271 x 600 MHz, X22 = 27r x 18 MHz, J = 271 x 200 MHz and the decay 
rate of the resonator is given by /c = 27r x 10 kHz. F is defined as the overlap between the ideally 
evolved density matrix and the digitally evolved one, F{t) = Tr{pi{t)pT{t)). The inset shows the 
overlap between the ideally evolved density matrix and the state of the system ait = 0,1 I\f2{c^ + 
! \/ 2 ) | 0 )p ® 111020304 )^, i.e., the same probability for having 1 and 2 photons in the resonator 
and all the superconducting qubits in the ground state of a? except the first, which is excited. 


Fig. 5 Protocol for the digital 
quantum simulation of a spin- 
chain coupled to a bosonic 
mode with superconducting 
circuits, in terms of unitary 
evolutions of Hamiltonians 
H\, H 2 and H, being the 
interaction defined in Eq. 
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3.3 Collective spin coupling mediated by resonators 

In this subsection, we extend the Ising model presented in Eq. <0 by adding three- 
body interactions. The method can be generalized to arbitrary interaction orders. 
This extension allows us to simulate problems such as quantum chemistry |[20l [36l 
[571 [38l, as well as fermionic lattice models EH EH, by using the Jordan-Wigner 
mapping to map fermionic interactions into spin interactions. The Hamiltonian for 
N spins including three-body interactions can be written as 

7/ = 7£cT?cr| + G^ cT?cr|cT/ + 5^cr^, (11) 

{jk) {jkl) j 

Here, we have added one collective interaction term with coupling constant G, 
which is the coupling among three nearest neighbour spins. This model can be 
simulated by enriching with additional gates the protocol for the Ising model in 
section |3.1| That is, together with single-qubit rotations along the v direction, 
Xy(0) = exp(—/0cTy), and two-qubit zz gates, ZZjk{6) = exp(—/0(Ty(7|), we also 
consider the combination of collective gates shown in Fig.|^ This will allow us to 
couple three qubits, ZZZjj^i{^) = exp(—/j8(7yCj|(7|). 

The collective spin interaction of this model can be decomposed into (A^ — 1) 
two-qubit gates and 2{N — 2) three-qubit gates, and the transverse field is mimicked 
by N single qubit rotations. Moreover, we notice that the digital error of the Trotter 
expansion in Eq. O is reduced due to the fact that the interaction summands of the 
Hamiltonian commute with each other. The Trotter expansion for this model reads 

e-iHt ^ t/s JZ{jic) I/s GZ{jki) t/s ^ ^j2) 


where 


exp(-/e £ a]al) = 

Uk) 

exp(—/j3 £ o^olof) = ... e~'^^N- 2 ^N-i^N^ 

Uki) ' 

cxp{-i(l) £ a;) = (13) 

j 

with 6 = —Jt, P = Gt and (j) = Bt, t being the simulated execution time. The col¬ 
lective gate for three qubits can be decomposed into two-qubit gates, as in Fig. 
Recently, the implementation of collective gates with a quantum bus has been pro¬ 
posed in superconducting circuits 15^ . 

In Fig.|7j we plot a numerical simulation of the extended Ising model with higher- 
order terms and transverse field for several Trotter steps. The figure shows as in the 
previous examples how the simulated dynamics with digital methods becomes more 
accurate when compared with the exact one when the number of Trotter steps is 
increased. 
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Fig. 6 a Protocol for performing one of the three-qubit interactions appearing in Eq. Ill} 
with collective gates. Here ZZZi 23 (a) = exp(—/aa^alal), = exp(—/0a^/2) is the rota¬ 
tion along the 7-axis of a qubit, X = exp(/aa^), and = exp(—a?aj). b The 

same interaction ZZZ can be realized with two-qubit gates where ZZa = exp(/7ra^ C) C7^/4), and 
ZZb = exp(—/Tra^ O a^/4). 



0 


Fig. 7 Digital simulation of the extended Ising model with a transverse field and three-body terms 
for four spins in superconducting circuits, where J = G = 2nx 400 MHz and B = 2nx 200 MHz, 
for a phase of 0 = —Jt = 4. The plot shows the fidelity of the digitally evolved state with the ideally 
evolved one for different Trotter steps, .y = 7,9,11. The inset shows the overlap between the ideally 
evolved state with the initial state, that is, all qubits in |0)^ state. 
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4 Conclusions 


In this article, we have proposed a digital quantum simulation of spin chains coupled 
to bosonic modes by means of circuit quantum electrodynamics architectures. We 
have presented a method for decomposing spin interactions and implementing them 
stroboscopically with available single and two-qubit gates. Furthermore, we have 
considered both circuit QED setups implementing capacitive couplings between su¬ 
perconducting qubits and transmission line resonators acting as quantum buses. We 
have exemplified our method with the quantum simulation of the Ising model with 
transverse field, a spin chain coupled to a bosonic field mode, and a many-body spin 
model with three-body terms, which are realized through a bosonic quantum bus. 
These results show that spin chains and bosonic field modes can be implemented 
efficiently with superconducting qubits. 

We thank Rami Barends for useful discussions and acknowledge support from 
Spanish MINECO FIS2012-36673-C03-02; Ramon y Cajal Grant RYC-2012-11391; 
UPV/EHU Project EHUA14/04, and two UPV/EHU PhD grants, Basque Govern¬ 
ment IT472-10; PROMISCE, and SCALEQIT EU projects. 
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